!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
! Initialisation parameters for shocktube problem
!
&init_pars
  cvsid='$Id: start.in,v 1.1 2013-06-17 12:04:56 wlyra Exp $',
  ip=14,                        ! debugging parameter
  xyz0  =  0.4,-3.141592,-0.048,       ! first corner of box
  xyz1  =  2.5, 3.141592, 0.048
!
  lperi = F  , T  , T,         ! periodic direction?!
!  bcx = 'a','a2','s','s','s','s' !'a2','s'
  bcx = 'out','a2','s','a2','s' !,'s' !'a2','s'
!  bcx = 'a','a2','s','a2','s'!,'s' !'a2','s'
!
  lwrite_ic=T
  coord_system='cylindrical'
  lcylindrical_gravity=T
!
  lwrite_aux=T
  llocal_iso=T
/
&initial_condition_pars
  density_power_law=1.5
  temperature_power_law=0.0
/
&eos_init_pars
  cs0  = 0.3162278
  rho0 = 1.
  gamma= 1.0
/
&hydro_init_pars
  inituu='gaussian-noise'
  ampluu=1e-3
/
&density_init_pars
  ldensity_nolog=T
  !initlnrho='gaussian-noise'
  !ampllnrho=1e-4
/
&grav_init_pars
  ipotential='no-smooth'
  g0=1.
/
&special_init_pars
/
&particles_init_pars
  initxxp='random-cylindrical', initvvp='Keplerian'
  dustdensity_powerlaw=1.5
  bcpx='rmv',bcpy='p',bcpz='p'
  lcommunicate_rhop=T
  lparticlemesh_tsc=T
  tausp=1.0
  eps_dtog=1.
  rp_int=0.6
  rp_ext=2.4
/
